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STABILITY  AND  SYMMETRY  IN  INERTIAL  CONFINEMENT  FUSION 


I.  Symmetry 

The  critical  consideration  for  achieving  high  spherical  compression  is 
to  determine  the  symmetry  requirements  on  the  driver  inposed  by  the  symmetry 
that  can  be  maintained  during  the  implosion  process.  Ablation  pressure  sym¬ 
metry  with  Ap/p  <  few  percent  oust  be  maintained  if  high  gain  is  to  be 
achieved1.  We  have  investigated  the  effect  of  an  asymmetric  laser  beam 
(1.05  pm)  with  an  intensity  ratio  Imax/Imin  ■  2  on  thin  foils 
(15  mn  CH)  with  our  FAST2D  code2.  We  find  that  lateral  energy  flow  in  the 
form  of  transverse  thermal  conduction  smooths  the  effects  of  laser  asymme¬ 
tries.  The  degree  of  smoothing  is  directly  related  to  the  distance  between 
the  ablation  and  critical  surfaces.  To  smooth  asymmetries  to  the  extent  that 
the  ablation  pressure  variation  is  less  than  a  few  percent  requires  that  the 
scale  length  of  the  asymmetry  be  comparable  to  the  distance  between  the  abla¬ 
tion  and  critical  surface. 

The  numerical  results  are  summarized  in  Fig.  1. ,  where  the  variation  in 
ablation  pressure  is  plotted  as  a  function  of  the  average  laser  intensities 
for  a  series  of  asymmetry  scalelengths  (A^)  in  the  range  of  100  pm  <  X^  < 

600  pm.  (Also  shown  is  the  equivalent  distance  between  the  ablation  surface 
and  the  critical  surface. )  These  results  are  in  agreement  with  the  NRL 
experimental  results3  shown  in  Fig.  2.  We  cotqpare  the  foil  velocity  nonuni¬ 
formity  -  1)  for  the  100  and  200  m  asymmetries  to  the 

experimental  results  for  a  1U0  pm  asymmetry  with  I^x^min  ■  2. 

A  new  analytic  steady-state  planar  ablative  flow  model **  has  been  devel¬ 
oped  from  which  scaling  laws  can  be  derived  which  depend  only  on  the  mate¬ 
rial,  laser  wavelength  and  absorbed  laser  irradiance.  The  derived  scaling 
Menuaeript  approved  September  20, 1982. 


1 


lavs  are  in  good  agreement  vlth  both  experiment  and  fluid  simulations.  On 
Fig.  2  we  show  the  results  of  the  analytic  model  for  laser  smoothing. 

IX.  Stability  -  Nonlinear 

The  Rayleigh~Taylor( RT)  instability  may  also  destroy  the  spherical  sym¬ 
metry  of  high  aspect  ratio  iiqploding  shells.  This  desymmetrizing  instability 
may  grow  in  the  ablation  layer  or  at  an  internal  interface  where  materials  of 
different  densities  abut,  thing  our  FAST2D  code  we  have  modelled  the  RT 
instability  far  into  the  nonlinear  regime  and  beyond  the  point  of  foil  frag¬ 
mentation  for  both  a  nultimode  perturbation5  and  for  3ingle  mode  perturba¬ 
tions  for  a  series  of  perturbation  wavelengths  (X)  in  the  range  of  1/2  <  X/AR 
<  10,  where  AR  is  the  cold  foil  thickness  (20  w).  laser  intensities  were  in 
the  range  l-5xio13  W/cm2. 

We  find  linear  growth  rates  well  below  classical  values  (by  a  factor  on 
the  order  of  3-1*)  and  a  cutoff  in  the  growth  rate  for  wavelengths  less  than 
the  cold  foil  thickness.  Hie  striking  result  is  the  dominance  of  nonlinear 
effects,  i.e. ,  the  Kelvin-Helnholtz  (KH)  roll-up,  for  short  wavelength  per¬ 
turbations.  Although  the  linear  growth  rates  increase  as  kl/2  up  to  the 
cutoff,  the  KH  roll-up  dominates  at  large  k,  drastically  reducing  the  pene¬ 
tration  rate  of  the  dense  spike  below  its  free  fall  value  and  effectively 
doubling  the  aspect  ratio  of  the  foil6. 

Hie  growth  rate  of  the  spike  shifts  from  exponential  to  quadratic  (free- 
fall)  when  the  anplitude  of  the  disturbance  is  on  the  order  of  1/5  the  wave¬ 
length  {A  0(0*2  X)].  As  a  result,  the  anqslitude  of  the  large  wavelength 
modes  (X/AR  >  1.25)  goes  into  "free-fall"  at  a  later  time  and  with  an  effec¬ 
tive  "free-fall"  acceleration  conparable  to  the  foil  acceleration  (0(l.5*1015 
cm/s2)l.  Hie  exponential  grcwth  of  the  short  wavelength  modes 


(X/AiR  <  1.25)  saturates  earlier  in  time  (since  the  linear  growth  rate 
is  larger)  and  the  "free-fall"  acceleration  of  the  spike  is  reduced.  This 
reduction  is  a  direct  consequence  of  the  tip-widening  brought  on  by  the  KH 
roll-up.  The  frontal  area  of  the  spike  is  increased  thereby  increasing  its 
drag  and  reducing  its  rate  of  fall. 

These  results  have  serious  implications  with  respect  to  shell  integrity 
as  a  function  of  perturbation  wavelength.  This  is  evident  from  Pig.  3  where 
both  the  "free-fall"  acceleration  of  the  spike  and  the  final  aspect  ratio  of 
the  foil  are  plotted  as  a  function  of  perturbation  wavelength  (X/AR).  Note 
that  the  "free-fall"  acceleration  ("g")  of  the  spike  is  reduced  by  over  a 
factor  of  two  for  the  short  wavelengths.  A  combination  of  the  reduced  "free- 
fall"  and  the  earlier  exponential  growth  saturation  increases  the  aspect 
ratio  (R/AR)  of  the  short  wavelength  perturbation  by  *  factor  of  two  over  the 
long  wavelength  modes.  Thus  nonlinear  effects,  in  the  form  of  the  KH  roll¬ 
up,  dominate  at  large  k,  saturating  the  linear  k1/2  form  of  thiA  RT  instabil¬ 
ity,  reducing  the  rate  of  "fall"  of  the  spike  and  increasing  the  aspect  ratio 
of  the  foil  before  breakup. 

III.  Stability  -  Linear 

We  now  show  that  the  vorticlty,  which  plays  such  a  dominant  role  in  the 
nonlinear  evolution  of  the  RT  instability,  is  also  responsible  for  the 
reduced  linear  growth  rates.  In  particular,  it  is  the  shedding  of  the  vorti- 
city  generated  at  the  ablation  layer  that  gives  rise  to  the  low  values  of  the 
growth  rate7. 

figure  U  shows  the  vortex  shedding  process  from  the  ablation  layer  for 
the  20  in  perturbation  wavelength  when  the  system  is  in  the  linear  regime. 

The  finite-sited  vortices  are  generated  by  the  baroclinic  (non-collinear 


density  and  pressure  gradients)  nature  of  the  profiles.  Uie  vorticity 
generated  at  the  ablation  surface  is  then  convected  with  the  flow  and  shed 
into  the  blovoff.  This  shedding  process  is  very  similar  to  the  shedding  of 
the  vortex  structure  that  arises  from  the  flow  about  a  bluff  body  and  the 
measured  Strouhal  number  (dimensionless  shedding  frequency)  agrees  very  well 
with  this  theory.  The  shedding  process  reduces  the  classical  RT  growth  rate 
by  removing  some  of  the  vortex  motion  from  the  growth  region.  We  express 
this  reduced  growth  rate  as  y*  *  (kg)1/2-af,  where  f  is  the  measured  shedding 
frequency  and  a  is  some  fraction  (less  than  one).  Figure  5  compares  this 
result  with  the  coiqputational  growth  rates. 

Not  only  does  this  shedding  reduce  the  linear  RT  growth  rate  but  a  sym¬ 
metric  von  Karman  vortex  street  develops  in  the  blow  off  region.  The  sym¬ 
metric  vortex  street  is  very  unstable  and  this  structure  could  have  serious 
effects  on  laser  absorption,  thermal  conduction  and  plasma  instabilities  in 
the  underdense  region. 

Vorticity  continues  to  be  shed  from  the  ablation  layer  until  the  spike 
amplitude  becomes  conparable  to  half  the  perturbation  wavelength  and  the 
bubble  region  ceases  to  ablate.  Hie  vorticity  is  then  advected  down  the 
sides  of  the  spikes  where  it  sets  up  a  reversed  flow  which  causes  the  tips  of 
the  spike  to  widen  thus  reducing  their  penetration  rate.  Vortex  dynamics 
plays  a  dominant  role  in  both  the  linear  and  nonlinear  phases  of  the  RT 
instability. 

We  nave  also  simulated  the  linear  and  nonlinear  evolution  of  the  RT 
Instability  for  strongly  perturbed  targets  and  these  results  are  presented 
elsewhere  at  this  conference6. 
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IV.  Three-Layer  RT  Problem 

The  interior-layer  Rayleigh-Taylor  problem  can  occur  at  many  different 
times  and  locations  during  the  iaplosion  of  a  complicated  nultiple-shell 
target.  The  preignition  mix  problem  is  the  most  familiar  situation.  This 
occurs  near  the  final  stages  of  cojipression  of  the  fuel  when  the  interface 
between  the  higher  density  pusher  and  the  fuel  compressionally  decelerates. 

A  three-layer  mix  problem  arises  where  a  light  fuel  layer  is  separated  from  a 
denser  pusher  shell  by  a  thin  but  very  dense  heat-shield  layer.  Short-wave¬ 
length  perturbations  permit  the  upper  and  lower  surfaces  of  the  heat  shield 
to  move  independently,  so  the  RT  instability  can  result  even  though  overall 
stability  is  ejq>ected  at  long  wavelengths.  In  order  to  investigate  the  non¬ 
linear  regime  of  the  RT  instability  for  these  two-  and  three-layer  mix  prob¬ 
lems  we  have  developed  a  Iagrangian  model  (SPLISH)  which  bases  its  fluid 
representation  on  a  two-dimensional  dynamically  reconnectable  grid  of  tri¬ 
angles  9. 

In  the  three  layer  mix  problem  the  densities  of  the  layers  are  in  the 
ratio  0. 1:2:1  so  that  the  unstable  interface  has  an  Atwood  number  of  1/3-  At 
this  Atwood  number  the  evolution  of  the  unstable  interface  is  governed  by  the 
interaction  between  the  RT  and  KH  instabiltie3  in  both  the  linear  and  non¬ 
linear  regimes.  The  evolution  of  the  RT  and  KH  instabilities  can  be  sepa¬ 
rated  into  four  distinct  phases:  (l)  the  initial  classical,  linear  growth 
phase,  ending  when  the  spike  to  bubble  height  is  half  the  perturbation  wave¬ 
length,  or  for  thin  layers,  when  the  spike  to  bubble  height  equals  the  layer 
thickness;  (2)  a  saturated  RT  growth  phase  with  steadily  increasing  KH  roll¬ 
up  and  Increasing  wavelength,  ending  when  the  spike  tip  becomes  wider  than 
the  spike  body;  (3)  a  snail  scale  structure  phase  and  (4)  a  turbulent  mixing 
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phase.  The  third  phase  is  characterized  by  downward  flow  from  the  thinned 
layers  above  the  bubbles  resembling  Jets  feeding  the  spike.  A  stagnation 
point  forms  at  the  point  where  the  Jets  collide,  forcing  an  upward  Jet 
opposed  to  the  spike  fall.  In  the  last  phase,  large  anplitude  waves  run 
across  the  stable  interface  as  well,  giving  rise  to  Jetting  and  wave¬ 
breaking.  Many  similarities  exist  between  the  ablative  and  nonablative 
cases.  Wie  KH  instability  plays  a  very  major  role  in  both  cases  in  that  it 
reduces  the  penetration  rate  of  the  dense  spike  below  its  free-fall  value. 

V.  Wavelength-Intensity  Window 

The  stability  and  symmetry  results  provide  limitations  on  the  designs  of 
the  laser  fusion  pellets  if  high  gain  is  to  be  achieved.  RT  results  indicate 
a  limit  in  the  aspect  ratio  of  around  ho.  Symmetry  results  indicate  that  the 
separation  distance  of  the  critical  from  the  ablation  surface  should  be  about 
.3  to  .5  of  the  current  pellet  radius.  We  have  used  our  one-dimensional  code 
(FAST1D)  to  investigate  scaling  laws  and  find  an  intensity  wavelength  window 
based  on  these  limits  10. 

Calculations  at  laser  wavelengths  of  0.53,  1.06,  and  2.7  l®  with  pellet 
radii  of  1,  2,  and  5  mm  were  used  to  derive  scaling  laws  for  the  pressure  P 
and  the  separation  ratio  0/R  versus  laser  intensity,  laser  wavelength,  and 
pellet  radius.  We  find  a  pressure  scaling  of  the  form 
P  ~  (IQ  b*25R-0«l8)0»7>  an<j  a  separation  ratio  scaling  of  the  form 
D/R  -  (laA3*8)°,7/(R+B).  In  the  region  of  interest,  B  *  1  mm,  and  is 
only  weakly  dependent  on  Ia  and  X.  This  form  goes  to  the  correct  limits 
for  both  large  and  small  radii. 

Our  symmetiy  results  indicate  that  if  the  separation  distance  D  is  more 
than  0.3  times  the  laser  inhomogeneity  wavelength,  then  6P/P  <  .ldl/I 
(FJ  1),  providing  sufficient  for  10%  laser  nonuniformities. 
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For  today's  moderate  aspect-ratio  fillet,  with  RQ/AR0  =  8-20,  an  abla¬ 
tion  pressure  of  about  10-30  Mbar  is  required  to  drive  the  shell  inward  to 


fusion-level  velocities,  with  use  of  a  relatively  unshaped  laser  pulse. 

Figure  6  shows  the  minimum  smoothing  and  pressure  barriers  based  on  these 
limits.  The  optimal  laser  intensity  is  approximately  a  few  times  1011*  W/cm2 
and  at  approximately  one-micron  laser  wavelength.  This  should  suffice  for 
producing  both  sufficient  ablation  pressure  and  sufficient  two-dimensional 
smoothing.  (For  a  flux  limit  of  f  =  0.03  the  optimal  wavelength  would 
increase  to  2-3  urn. ) 

For  fixed  smoothing  decreasing  the  laser  wavelength  X  increases  the 
instability  parameter  IX2  increasing  the  potential  for  plasma  instabilities. 
Equally  important,  the  high  ablation  pressure  at  large  intensities  might  be 
too  large  thereby  driving  strong  shocks  that  would  preheat  the  pellet. 

There  3eems  to  be  a  barrier  at  shorter  wavelengths  because  the  smoothing 
requirement  quickly  drives  one  to  unacceptably  high  laser  intensities  and 
ablation  pressures.  And  at  longer  wavelengths,  the  ablation  pressure 
requirement  on  the  laser  intensity  soon  exceeds  the  threshold  for  plasma 
instabilities.  The  actual  location  of  the  wavelength  window  is  somewhat 
uncertain,  but  it  appears  to  lie  between  1/2  and  2.7  ym. 

Another  important  consideration  in  parameter  mapping  is  the  hydrodynamic 
efficiency.  In  Fig.  6  we  also  show  the  variation  of  intensity  with  laser 
wavelength  for  fixed  hydrodynamic  efficiency.  Efficiencies  greater  than  7% 
are  achieved  below  this  line.  The  highest  efficiency  is  achieved  by 
operating  at  the  intersection  of  the  smoothing  requirement  and  the  pressure 
limitations. 


Fig.  3  —  A  plot  of  the  “free-fall"  acceleration  (“g”  in  units  of  1014  cm/s2)  of  the 
spike  and  the  final  aspect  ratio  of  the  foil  as  a  function  of  the  perturbation  wavelength 
(X/AR).  The  solid  circles  are  the  aspect  ratios  obtained  with  a  laser  intensity  of  5  X 
1013  W/cm2  and  the  #  symbol  is  the  aspect  ratio  obtained  with  a  perturbation  wave¬ 
length  of  200  Mm. 
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Fig.  4  —  Vorticity  contours  for  the  20  pm  mode  (X/AR  *  1)  when  the  system  is  in  the 
linear  regime.  The  contours  are  in  10%  increments  of  the  maximum  vorticity  (+2.75  X 
108,  —  2.87  X  108  S'1).  The  solid  (dashed)  contours  are  positive  (negative)  vorticity. 


The  heavy  lines  delineate  the  10%  (right)  and  90%  (left)  density  contour. 
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Fig.  5  —  A  comparison  of  the  computational  growth  rates  (— < e— )  to  the  classical  value 
((kg)**. — )  and  to  the  vortex  riwdding  theory  (0).  (a  -  0.56  to  fit  the  20  »m  mode). 
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LASER  WAVELENGTH  (#<m) 

Fig.  6  —  Wavelength-intensity  window  determined  by  the  symmetry  requirements 
(downward  sloping  solid  curves)  and  the  stability  limits  (aspect  ratio  <  40)  for  two 
radii  pellets  (.1  and  .5  cm).  Constant  efficiency  of  7%  is  shown  by  dashed  line.  Opera¬ 
ting  region  is  the  upper  right  hand  quadrant. 
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